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The NA49 fixed-target experiment studied high energy-density matter pro- 
duced in nucleus-nucleus reactions at the CERN SPS. In central Pb-nPb colli- 
sions at 158yl GeV the energy density at the early stage substantially exceeds the 
threshold for quark deconfinement predicted by lattice QCD. The produced mat- 
ter shows strong transverse and longitudinal flow. Ratios of yields of produced 
particles are approximately consistent with statistical equilibration. An energy 
scan through the SPS range revealed structure in the energy dependence of tt and 
K yields as well as of the inverse slopes of transverse mass distributions. These 
features suggest that a deconfined phase starts to be produced at around 30 A GeV 
in central Pb-nPb collisions. The analysis of fluctuations and correlations has not 
yet provided evidence for the predicted critical point of QCD. 
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1. Introduction 

Qualitative considerations based on the finite size of hadrons IT] as well as 
quantum chromodynamics (QCD) calculations on the lattice f2l f3l predict that at 
sufficiently high energy density strongly interacting matter will transform into a 
state of quasi-free quarks and gluons, the quark gluon plasma (QGP). Theoretical 
investigations also found |4] that this phase transition is of first order for finite 
quark masses and large non-zero baryon density. The phase boundary is predicted 
to end in a critical point and turn into a rapid crossover as the net baryon density 
decreases. 




Fig. 1. Schematic layout of the NA49 experiment at the CERN SPS showing beam detec- 
tors, superconducting dipole magnets, time projection chambers (VTPC, MTPC), time-of- 
flight arrays (TOP) and calorimeters (RCAL, VCAL). A thin solid target T is used for A+A 
collisions (a), which is surrounded by a detector of slow protons (CD) for p+A collisions 
(c). A liquid H2 target is employed for p+p collisions (b) 

The initial stage of high energy collisions of large nuclei provides the best en- 
vironment to produce the deconfined phase of matter in the laboratory [5,, 6]. Lead 
ions first became available at the CERN SPS in 1994. In the first publication from 
this programme [7| NA49 demonstrated that in central Pb-i-Pb collisions at top 
SPS energy the initial energy density exceeds the critical value of « 1 GeV/fm^. 
The SPS experiments found that the reactions produced an explosively expanding 
fireball. Moreover, originally proposed signatures of the QGP, i.e. J/^* suppres- 
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sion, strangeness enhancement, and possibly thermal photons and dileptons were 
observed [8]. However, these signatures are not specific for deconfinement. The 
NA49 collaboration therefore performed an energy scan from 20 - 158^ GeV in 
order to search for structure in the energy dependence of hadron production charac- 
teristics which could indicate the onset of deconfinement [9]. The measurements 
indeed suggest structure around 30^4 GeV ifTOl which will be discussed below. 
NA49 also searched for fluctuations which might occur if distinct phases coex- 
isted in the early stage of the reactions or if hadrons froze out close to the critical 
point estimated LIJ to lie in a region accessible at the SPS. 

2. Experiment NA49 at the CERN SPS 

The main features of the NA49 experiment [11] located in the H2 beam line 
of the North Experimental Hall (see fig. [Hi are large acceptance precision tracking 
(Ap/p^ (0.3 — 7) • 10~^(GeV/c)~^) and particle identification in the central and 
forwards rapidity regions using time projection chambers (TPCs). Charged parti- 
cles (vr, K, p, p) are identified mostly from the measurement of their energy loss in 
the TPC gas (accuracy 3-6 %). Yields are obtained by fitting a sum of Gaussian 
functions for the various particle species to the dE/dx distributions in small bins of 
momentum p and transverse momentum px- At central rapidity the identification 
is further improved by measurement of the time-of-flight (resolution 60 ps) to ar- 
rays of scintillation counter tiles. Strange particles (K^, A, E, Q) are detected via 
decay topology and invariant mass measurement. A forward calorimeter measures 
the energy of the projectile spectators from which the impact parameter in A-i-A 
collisions and the number of participating nucleons are deduced. Also reactions 
of C and Si ions as well as of deuterons were studied. These beams were pro- 
duced by fragmentation of the primary Pb beam and selected by magnetic rigidity 
and by specific energy loss in transmission detectors in the beam line. In addition 
NA49 pursued an extensive program of proton induced collisions for a study of the 
evolution of particle production from p-i-p via p-i-Pb to Pb-i-Pb reactions. 

3. Analysis procedure 

All events used for physics analysis were required to have a good vertex in the 
respective target. Tracks of accepted particles had to fulfill quality requirements 
on the number of measured points and the distance of closest approach to the 
event vertex. Particle yields have been corrected for acceptance, reconstruction 
efficiency and feed-down contamination from weak decays. Results on particle 
yields for the upper 3 energies have been published in |[T2l [T3l [141 where more 
details on the analysis procedures may be found. Results at 20 A and 30 A GeV are 
preliminary. 
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Fig. 2. Transverse mass mj spectra at midrapidity in central Pb+Pb collisions at 20A 
(left), 40A (center) and 158^ (right) GeV. The curves show the result of a blast wave fit 
lfT6l with parameters: temperature If and surface radial flow velocity po- 



4. Thermal freezeout parameters from spectra and correlations 

Midrapidity invariant yields as a function of the transverse mass mt are shown 
in fig. |2] at 20 A, 40 A and 158^4 GeV as examples for the large variety of par- 
ticle species measured by NA49. The spectra become progressively flatter with 
increasing particle mass, a fact that can be explained by the combined effect of 
the random thermal momentum distribution and strong radial flow in the produced 
matter droplet. The development of collective flow is generally attributed to the hy- 
drodynamic pressure generated in the dense early stage [ 15|. A hydrodynamically 
based "blast wave" parameterization fT6l indeed provides a reasonable description 
of all the spectra with two parameters: a temperature w (90 - 1 10) MeV and a 
radial flow velocity of /Oq ^ 0.8 c at the surface. These parameters characterize the 
thermal/kinetic freezeout of the fireball. The pion spectrum overshoots the model 
curve at low mt due to the feed-down contribution from resonance decays and this 
region was therefore not used in the fit. More detailed analyses ifTTl [141 seem to 
suggest an earlier freezeout of multistrange hyperons at higher temperature and 
smaller radial flow velocity. 

Collective flow also affects quantum statistics induced correlations of identical 
particles. For negative pion pairs the range of these Bose-Einstein (BE) corre- 
lations in momentum difference, which is inversely proportional to the effective 
pion source size, is reduced by flow. This reduction gets stronger with increasing 
transverse momentum (see fig. [3] (left)). In the NA49 publication |18| this effect 
together with the nit spectrum of pions was used for the first time to demonstrate 
strong collective radial flow in Pb-nPb collisions, i.e. the correlation of the radial 
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Fig. 3. Left; Projections of the tt^tt^ correlation function on momentum difference q 
in the transverse (side, out) and beam (long) directions for central Pb+Pb collisions at 
15SA GeV. Right: Effective radius parameters R of the pion source versus mean transverse 
momentum kt of the pion pairs in central Pb+Pb collisions at the five SPS energies. The 
curves show results of simultaneous fits of the blast wave model to the radius parameters 
and the transverse mass mt spectra of pions and protons. All results refer to midrapidity 
pion pairs. 



position of the emission points with the transverse momentum of the produced 
particle. This information cannot be obtained from transverse mass spectra alone. 
The fitted effective source sizes, the so-called radius parameters R, are shown as 
a function of averaged transverse momentum kt of the pion pairs for all the SPS 
energies in fig. [3] (right). The decrease of R with increasing kt due to radial flow 
is clearly observed. On the other hand, as shown in fig. |4] (left) the radius parame- 
ters show remarkably little energy dependence over the whole energy range from 
AGS to RHIC. In particular, there is no indication of an increase of Rout at the 
SPS which might be expected from a first order phase transition. In fact, hydro- 
dynamic calculations with a simple freezeout procedure overpredict the values of 
Rout and Riong (15). A solution of this so-called HBT puzzle seems to lie in a more 
sophisticated treatment of the frezzeout process [,19] . 

The hydrodynamics inspired blast wave parameterisation provides a compre- 
hensive description of the kinetic freezeout stage of nucleus-nucleus reactions. 
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The parameterization fT6l assumes a uniform pion emission density in a cylinder 
of radius R and a radial flow velocity increasing linearly with radius to a surface 
maximum of p^. Further parameters of the model are the kinetic freezeout tem- 
perature T (mainly determined by the rrit spectra), the lifetime r of the fireball 
(derived mostly from Riong) and the duration At of the pion emission burst (ob- 
tained from the difference of Rout and Rside)- Resulting parameter values from 
simultaneous fits to radii and mt spectra from AGS to RHIC energies are plotted 
versus collision energy in fig.|4](right). The corresponding radius parameters at the 
SPS are shown as curves in fig. [3] (right). The fit results demonstrate that the mat- 
ter droplet at freezeout has increased radially by a factor of two and has attained 
strong collective radial flow. While frezeout temperature and fireball lifetime ap- 
pear to to increase slowly with energy, the surface radial flow velocity seems to 
saturate at SPS energies. 








ft.*" 








AGS (E89S) 
SPS (NA49) 
RHIC (STAR) 



10 10- 

Fig. 4. Left; Energy dependence of radius parameters R of the pion source at mean trans- 
verse momentum kt= 0.2 GeV/c and midrapidity in central Pb+Pb and Au+Au collisions. 
Results from AGS, SPS and RHIC are shown as squares, filled (NA49) and open (NA45) 
circles, and triangles respectively. Right: Freezeout parameters obtained from simulta- 
neous fits of the blast wave model |fT6l to radius parameters from tt^tt^ Bose-Einstein 
correlations and mt spectra of pions and protons at midrapidity in central Pb+Pb (Au+Au) 
collisions from AGS (blue dots) through SPS (red dots, NA49 data) to RHIC (green dots) 
energies. 



While radial flow manifest in mt spectra and BE correlations builds up over 
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the full lifetime of the fireball, anisotropic flow generated in non-central collisions 
is particularly sensitive to the properties of the produced matter during the early 
phase of the reaction. In the almond shaped interaction region anisotropic pressure 
gradients and rescattering lead to an azimuthally anisotropic momentum distribu- 
tion. Since the spatial anisotropy disappears quickly, the momentum anisotropy 
is generated predominantly at early times in the reaction. The pressure rise with 
initial energy density is qualitatively expected to slow down when deconfinement 
sets in. Anisotropic flow is quantified by the Fourier expansion coefficients of the 
azimuthal angular distribution with respect to the reaction plane. At midrapidity 
the dominant effect is observed in the second coefficient V2, termed elliptic flow. 




Fig. 5. Left: Azimuthal distribution of A hyperons with respect to the reaction plane as 
determined from pions. Full dots show data, open dots mixed events. Right: V2 of tt, p, A 
versus transverse momentum at midrapidity in Pb+Pb collisions at 158A GeV. Curves 
show results from the blast wave model lfT6l . 

NA49 measured anisotropic flow of pions and protons in Pb-i-Pb collisions at 
40^ and 158 A GeV [20] and of A hyperons ^ at 158 A GeV. An example of 
the A azimuthal distribution with respect to the reaction plane as determined from 
pions is displayed in fig. [5] (left). The elliptic flow coefficient V2 at midrapid- 
ity is plotted versus transverse momentum px for the different particle species in 
fig- 13 (right). One observes a strong increase with pT and a clear mass hierarchy 
as predicted by hydrodynamic models fT5l. Values of V2 are about 30 % higher 
at RHIC (not shown) than at SPS. Both at the SPS and at RHIC the observed el- 
liptic flow V2 can be reproduced by the blast wave model parameterisation fl6l 
with parameters consistent with those describing spectra and Bose-Einstein cor- 
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relations f2Tl. Pure hydrodynamic model calculations overpredict V2 particularly 
at SPS energies. However, coupling of a hadronic rescattering phase to the initial 
deconfined hydrodynamic phase results in reasonable description of the measured 
energy dependence from SPS to RHIC |[22l . 

20AGeV 30AGeV 40AGeV 80 AGeV 158AGeV 
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Fig. 6. Rapidity distributions in central Pb+Pb collisions at 20^1 - 158 A GeV. SoUd dots 
show measurements, open dots were obtained by reflection around midrapidity. 



5. Hadron freezeout parameters from particle yields 

The large acceptance of the NA49 detector allows measurements of rapidity 
spectra from midrapidity up to almost beam rapidity (see fig.©. Rapidity distribu- 
tions of most particles are peaked at midrapidity. Only A hyperons, sharing a va- 
lence quark with the projectile, show a flattening at higher collision energy. Due to 
the reflection symmetry of Pb-i-Pb collisions 4ir yields can be determined. Particle 
yield ratios in A-i-A (as well as elementary) collisions are consistent with statistical 
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model predictions from threshold to the highest energies using only 3 parameters: 
a temperature T, a baryo-chemical potential fi b and a strangeness suppression pa- 
rameter These parameters characterize the freezeout of particle composition 
after which only the momentum distributions further evolve via elastic scattering. 
An example of the results of the fit with the statistical model for central Pb + Pb 
collisions at 158^ GeV US is shown in fig. H (left). 




Fig. 7. Left: Results of a statistical model fit to NA49 hadron yields in 158A GeV central 
Pb + Pb collisions [23 1. Right: Phase diagram of hadronic matter showing the hadron 
composition freezeout points obtained from statistical model fits ||23] . The solid curve 
represents the empirical freezeout condition 1 GeV/hadron [24]. The grey band indicates 
the expected first order phase boundary between QGP and hadrons which ends in a critical 
point E estimated from lattice QCD (|3l- 

The resulting parameters T and fiB f23l| are plotted in the phase diagram of 
hadronic matter in fig. |7] (right) together with the phase boundary predicted by 
lattice QCD [3|. One observes that the freezeout points at SPS energies approach 
the phase boundary and and are closest to the critical point at about 60^1 GeV. 



6. Indications for the onset of deconfinement 

A detailed overview of the energy dependence of strangeness production is 
presented in fig. [H The (K+)/(7r^") ratio (left) shows a steep rise from the thresh- 
old of kaon production, a maximum around the lowest SPS energy and a decrease 
to a somewhat lower plateau value. Although the microscopic transport models 
RQMD, UrQMD |25 1 and the statistical hadron gas model (full equilibrium ver- 
sion (7s = 1) supplemented by the freezeout condition 1 GeV/hadron l24l ) follow 
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\Is (GeV) ^ (GeV) ^ (GeV) 

Fig. 8. Ratio of Air yields versus nucleon-nucleon cms energy y^: (K^) / (tt^) (left); 
{K~)/{-K~) (center); (A)/(0.5 • ((7r+) + (tt"))) (right). Results from p+p collisions 
(open dots) and model predictions for Pb+Pb collisions (curves) are also shown. 

the gross trend, they do not reproduce the pronounced peak of the (K+)/(7r+) ratio 
at the SPS. Since anti-baryon production yields are small, (K+) counts essentially 
half of all s quarks in the final state (the other half is contained in K'^). In con- 
trast s quarks are distributed between anti-kaons and hyperons (mainly A) because 
of the large net baryon density at SPS energies. As a consequence of the rapid 
decrease of the net baryon density over the SPS energy range, the sharp peak in 
(K+)/(7r+) is reflected in a break in the energy dependence of the {K^)/{tt^) ratio 
(fig. [8](center)) and a wider maximum in the (A) / (tt) ratio (fig. |8](right)). 




Fig. 9. Left: Estimated multiplicities of all s and s valence quarks contained in produced 
particles as a function of energy. The straight line was obtained from a fit to the low 
energy points. Right: Ratio of this quantity to the pion multiplicity. The full curve shows 
the prediction of the statistical model using a parameterization of /iB(^/s) and 7s = 1 
ll24l . The dashed curve depicts the prediction of the SMES model (91 . 

Using isospin symmetry as well as predictions from the hadron gas model 
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for some unmeasured strange particle species one can obtain an estimate of the 
number of s and s valence quarks contained in the produced particles. The re- 
sult shown in fig. |9] (left) demonstrates that there is a break in the energy depen- 
dence of total strangeness production at 30 A GeV. Dividing by the pion multiplic- 
ity (fig. |9] (right)) confirms the sharp peak already seen in the (K+)/(7r^) ratio. It 
has been pointed out that these features coincide with a transition from a baryon 
to a meson dominated final state |26|. However, as seen from the statistical model 
curve in fig. |9] (right) this effect cannot adequately describe the measurements. 




Fig. 10. Left: Relative strangeness yields E5 — ((K + K) + {A))/{tt) versus Fermi 
energy variable Fw s^tv- Center: Pion yield per participating nucleon {Tr)/{Nw) versus 
F. Right: Energy dependence of the effective temperature parameter T of K+ transverse 
mass spectra. Results from p+p collisions are shown as open circles, model predictions 
for Pb+Pb collisions as curves. 

A close approximation of the total strangeness to pion production ratio can also 
be calculated from measured yields via the observable E5 = ((K + K) + (A))/ (vr) 
which is plotted in fig. [TO] (left) as a function of collision energy. As expected, it 
shows the same sharp peak followed by a plateau as the (K+)/(7r+) ratio, which 
is not seen in p-i-p collisions and not reproduced by hadronic models. On the 
other hand this feature can be understood in a reaction scenario with the onset of 
deconfinement around 30^ GeV as proposed in the statistical model of the early 
stage (SMES [9], dash-dotted curve in fig. [TOl (left)). reflects the behavior of 
the ratio of the number of s -1- s quarks to entropy in the model. It shows a steep 
threshold rise while the system stays in the hadron phase but drops to the value 
expected in the QGP at higher collision energies when the fireball initially reaches 
the deconfined partonic phase. At the same transition energy of 30^ GeV the rate 
of increase of the number of produced pions per participating nucleon (a measure 
of the entropy per baryon in the model) was predicted to increase. As seen in 
fig- [El (center), this seems to be confirmed by the measurements for collisions of 
heavy nuclei, while there is no change in energy dependence for p-i-p reactions. 

The transverse mass spectrum of kaons is well described by an exponential 
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Fig. 11. Energy dependence of {nit) — mg for 7r~ and tt"*" (left), K+ and K~ (center), and 
p and p (right) from AGS, SPS and RHIC experiments. Full (open) symbols show results 
for negatively (positively) charged particles. 



function e~™'/-^. As shown in fig. [TO] (right) for K+ the inverse slope parameter T 
rises steeply at low energies, stays at an approximately constant value through the 
SPS energy range and ends up slighly higher at RHIC energy. Similar behavior 
is seen for K~ (not shown), but is not observed in p+p collisions. These features 
have been attributed [27 ] to the constant pressure and temperature when a mixed 
phase is present in the early stage of the reaction. In fact, a hydrodynamic calcu- 
lation |[28ll modeling both the deconfined and the hadronic phases can provide a 
quantitative description (dash-dotted curve in fig. [TO] (right)). The spectra of other 
abundantly produced particles are not well fitted by exponential functions (see e.g. 
fig. O but their shape can be characterised by the average value of the transverse 
mass {mt) minus the rest mass mo. The results in fig. [TT] show that the near con- 
stancy of the slope of the spectra in the SPS energy range is a common feature for 
all these particles species. 



7. Results on fluctuations 

If the fireball freezeout occurs close to the boundary of a first order phase 
transition or near the QCD critical point then large event-to-event fluctuations may 
be expected. The NA49 detector was therefore designed with a large acceptance 
to allow for a wide range of fluctuation and correlation measurements. 

The average transverse momentum (pt) in the event is related to the temper- 
ature and radial flow in the fireball. The latter might fluctuate with the initial 
energy density and the fraction of deconfined matter at the early stage of the col- 
lision. The distributions for event- wise (pt) are shown in fig. [12] (left) for central 
PbH-Pb , Si-i-Si and p-i-p coUisions at 158^4 GeV by the data points |29|. These 
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Fig. 12. Left: Distributions of event-wise {px) in the rapidity interval !.!<?;< 2.6 for 
Pb+Pb, Si+Si and p+p collisions at 158A GeV. Right: Distribution of the event-wise K/tt 
ratio in central Pb-nPb collisions at 158/1 GeV for laboratory momenta > 3 GeV/c. Results 
from real events are shown by data points from mixed events by histograms. 



are compared to results from mixed events (histograms) in which correlations are 
destroyed by construction while preserving the inclusive px and multiplicity distri- 
butions. Clearly, the fluctuations are close to the statistical limit as represented by 
the mixed events and no distinct event classes are observed. A quantitative analysis 
was performed using the quantity <l>pj, which is defined to vanish for purely statis- 
tical fluctuations and which is furthermore independent of the number of coUison 
participants in superposition models. The results for <I>py in the rapidity range 
1.1 < < 2.6 are plotted versus collision system size in fig. [13] (left). The 
measured <I>py values are smaller than 5 MeV/c and thus below a few percent of 
the (pt)- Interestingly a significant system size dependence is observed with a 
maximum for peripheral Pb-i-Pb collisions. A similar dependence was found for 
the fluctuations of the multiplicity of negatively charged particles [,10,| . The two 
observations may be related due to a possible correlation between (px) and multi- 
pUcity ffl. 

The energy dependence of the fluctuation measure <^p^ at the SPS for cen- 
tral Pb-i-Pb collisions is plotted in fig. [13] (right). The observed values are small 
and approximately constant. Somewhat larger results have been obtained near 
midrapidity by experiment NA45 I3TI . At RHIC the {pr) fluctuations increase 
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Fig. 13. Left; The {pt) fluctuation measure for negatively chai-ged paiticles versus 
the number of nucleon participants Nw in p+p, C+C, Si+Si and Pb+Pb collisions at 158 A 
GeV. The curve shows predictions from the HIJING model filtered through the acceptance 
of the NA49 detector. Right (preliminary): Energy dependence of {pr) fluctuation mea- 
sure for negatively charged particles in central Pb+Pb collisions. The curve shows 
predictions from the UrQMD model filtered through the acceptance of the NA49 detector. 
Results are for the rapidity interval l.\ < y-^ < 2.6. 



substantially due to the rise of jet production ll32l . 

Fluctuations of the K/vr ratio are believed to increase near the critical point. It 
has also been pointed out [33] that at the onset of deconfinement there could be a 
significant decrease. The majority of strange quarks in the produced particles are 
contained in kaons at SPS and RHIC energies. The mass of strangeness carriers 
in the hadron phase (kaons) is much larger than the temperature whereas in the 
deconfined phase (strange quarks) it is smaller. Therefore fluctuations in the initial 
energy density at fixed collision energy are expected to lead to large fluctuations 
of the K/vr ratio when the system stays in the hadronic phase. However, these K/tt 
ratio fluctuations should significantly decrease if deconfinement occurs during the 
early stage of fireball evolution |j9l[33l. The SMES model curve in fig. [TOl (left) 
provides an illustration of this generic prediction when one looks at the change of 
the K/vr ratio for small variations of the energy variable F in the respective regions. 

The {K^ + K^) / {t:^ +71^) (as well as {p+p) / {-k'^ +tt~)) ratio was extracted 
event by event with a maximum likelihood fit method for particles with laboratory 
momenta larger than 3 GeV/c in central Pb+Pb collisions at the five SPS energies 
ll34l . The momentum cut is necessary to allow identification by specific energy 
loss of the particles in the TPCs. This leads to acceptances which shift to larger 
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Fig. 14. Energy dependence of dynamical fluctuations of the {K^ + K )/(7r+ + tt ) 



ratio (left) and the (p 
(preliminary) 



- p) / (7r+ + TT ) ratio (right) in central Pb+Pb collisions at the SPS. 



cms rapidities with increasing particle mass. An example of the distribution of 
the + K")/(7r+ + vr") ratio at 158A GeV is shown in fig. [H] (right) where 
it is also compared to the results from mixed events (histogram) [35|. Clearly, 
both distributions are very similar and indicate that non-statistical fluctuations are 
small. Quantitatively, dynamical fluctuations were calculated by comparing the 
width (J^ata of the ratio distribution to the one obtained for mixed events from 
the formula = g\^^^ - cr^-^g^. As seen from fig. [14] (left) the dynamical 
fluctuations of the (K^ + i^~)/(7r+ + vr") ratio amount to only around 5 % of 
the average ratio (about 0.15), but they show a decrease towards higher energy. 
Preliminary results from the STAR experiment at RHIC for higher energies also 
indicate values close to those at the highest SPS energy [36,1 . Such a behaviour 
is not seen in typical reaction models, e.g. UrQMD, but is consistent with the 
qualitative expectations from the onset of deconfinement fST]. 

The dynamical fluctuations of the (p + p)/ (vr+ + vr^) ratio are plotted for the 
five SPS energies in fig. [14] (right). The observed negative values are described 
quantitatively by the UrQMD model and are most likely due to the correlated pro- 
duction of p and vr via baryon resonances. 

Electric charge fluctuations were predicted to be strongly reduced if the fireball 
passed through a QGP phase [37 1. Fig. [15] (left) depicts the distribution of net 
charge in the papidity range —1.4 < y < 2.4 for central Pb-i-Pb collisions at 
158 A GeV as an example. As a consequence of charge conservation the data (full 
histogram) display a narrower distribution than mixed events in which this effect 
is removed. Fig. [15] (right) shows the measurements of electric charge fluctuations 
in terms of the variable A<I>q [38 | in which effects of global charge conservation 
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Fig. 15. Left: Distibutions of net charge Q = N+ - for the papidity range —1.4 < 
y < 2.4 in central Pb+Pb collisions at GeV for data (full histogram) and mixed 

events (dotted histogram). Right: Measure A"I>q of electric charge fluctuations corrected 
for charge conservation versus rapidity acceptance. 



are subtracted. Neither significant energy dependence nor the predicted reduction 
are observed. However, more recent model calculations suggest that hadronisation 
effects f39l and at SPS also resonance decays l'40'l probably mask the reduction. 

A first order phase transition from the deconfined to the hadronic phase during 
the early stage of Pb-i-Pb collisions is expected to result in a softening of the equa- 
tion of state and consequently in a long lifetime of the mixed phase. The study 
of the balance function (BF) was suggested as a method to investigate the time of 
hadron formation BTIl . In many models pions are assumed to be created in charge 
neutral pairs. Due to the buildup of longitudinal flow and rescattering during the 
evolution of the fireball, the rapidity correlation of opposite charges will be diluted 
for pairs formed early in the process but will remain narrow for late hadronisation 
as expected in a mixed phase scenario. This correlation is quantified by the BF 
which is plotted as a function of the pseudorapidity difference Ar] of the pairs in 
fig.[T6](left). The data are compared to shuffled events, in which the pseudo rapidi- 
ties of the particles are scrambled while the charges are retained. In this way one 
obtains an estimate of the maximum decorrelation under the condition of charge 
conservation. It is evident that the data show increasingly shorter range opposite 
charge correlations beyond the effects of charge conservation as the centrality of 
Pb-(-Pb collisions increases. 

The widths {A-q) are calculated numerically from the BF distributions and the 
results are plotted in fig. [T6](right.top) for the top SPS energy |42| and for RHIC 
fl3l . The width decreases by 17 % respectively 14 % for central compared to pe- 
ripheral colUsions for data while it remains constant for shuffled events. Thus the 
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Fig. 16. Left: Balance function versus pseudo-rapidity difference of oppositely charged 
particle pairs at midrapidity and top SPS energy. Curves show Gaussian fits to data (solid), 
shuffled events (dotted) and HIJING simulation (dashed). Right, top: Width of the balance 
function versus impact parameter at midrapidity in Pb+Pb collisions from NA49 at SPS 
(dots) and STAR at RHIC (stars). Solid symbols denote data, open symbols show results 
from shuffled events. Right, bottom: Energy dependence of the normalized narrowing 
parameter W of the balance function. 

predicted narrowing is observed, although it is now believed that a large part of 
this effect is due to the increasing strength of radial flow. In fact a model assum- 
ing hadron formation by quark coalescence as well as collective radial flow can 
reproduce the observed width of the BF [441. 

In order to be able to compare measurements at different energies and with dif- 
ferent acceptance one may calculate the normalized difference (in percent) W = 100- 
(^^^^sh«//^ - {Ar])))/{Ar]'''''ff) of the widths of the BF for shuffled and real 
events. The energy dependence of the normalised narrowing parameter W is plot- 
ted in fig. [T6](right.bottom). It shows a steady rise, i.e. more narrowing of the BF 
with respect to uncorrelated particle production throughout the SPS energy range 
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up to RHIC. One does not observe a significant structure at 30 A GeV, wiiere otlier 
observables suggest tlie onset of deconfinement. 

8. Conclusions 

Results from a comprehensive study of nucleus-nucleus collisions in the SPS 
energy range were presented. The produced matter shows strong transverse and 
longitudinal flow. Ratios of yields of produced particles are approximately consis- 
tent with statistical equilibration. The freezeout parameters of the hadron compo- 
sition are close to the phase boundary predicted by QCD. The main charcteristics 
of nucleus-nucleus collisions at the highest SPS energy are quite similar to those 
found at RHIC. 

The energy scan of central Pb-i-Pb collisions revealed that at around 30^4 GeV 
the ratio of strangeness to pion production shows a sharp maximum, the rate of 
increase of the produced pion multiplicity per wounded nucleon increases and the 
effective temperature of pions and kaons levels to a constant value. These features 
are not reproduced by present hadronic models but find a natural explanation in a 
reaction scenario with the onset of deconfinement in the early stage of the reaction 
at low SPS energy. Further support for this scenario is provided by the decrease of 
fluctuations of the K/vr ratio in the SPS energy range. On the other hand, measure- 
ments of fluctuations of the average transverse momentum, of the electric charge, 
as well as of correlations such as the balance function and vrvr Bose-Einstein cor- 
relations exhibit a smooth energy dependence. 

The present results suggest that further study of the role of the initial volume in 
the onset of deconfinement and the search for the predicted critical point of QCD 
are of great interest. This led to an initiative for a low energy program at RHIC 
P31 and the proposal of a future light-ion program at the SPS P6ll47l . 
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